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Abstract 

We consider the Dirichlet boundary value problem for nonlinear sys- 
tems of partial differential equations with p-structure. We choose two rep- 
resentative cases: the "full gradient case", corresponding to a p-Laplacian, 
and the "symmetric gradient case", arising from mathematical physics. 
The domain is either the so called "cubic domain" or a bounded open sub- 
set of K 3 with a smooth boundary. We are interested in regularity results, 
up to the boundary, for the second order derivatives of the velocity field. 
Depending on the model and on the range of p, p < 2 or p > 2, we prove 
different regularity results. It is worth noting that in the full gradient 
case, with p < 2 , we cover the degenerate case, and obtain W^ 2,9 -global 
regularity results, for arbitrarily large values of q. 

Keywords: p-Laplacian systems, regularity up to the boundary, full regularity. 



1 Introduction 

We are concerned with the regularity problem for solutions of nonlinear systems 
of partial differential equations with p-structure, p > 1, under Dirichlet bound- 
ary conditions. In order to emphasize the main ideas we confine ourselves to 
the following representative cases (where \l > is a fixed constant): 
The "full gradient case" 

(1.1) -V-S(Vti)=/, 
where 

(1.2) fiT(Vtt) - (m+ I Vu\f- 2 V«; 
and the "symmetric gradient case" 

(1.3) -V-S(Vu)=f, 
where 

(1.4) S{Vu) = (fi+ \Vu\f- 2 Vu. 
As usual, 

Vu = i(Vn + V/) 

is the symmetric part of the gradient of u. 

When fi = in (|1.2[> , the system (jl.ll) is the well-known p-Laplacian system. 
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It is worth noting that our results concern global (up to the boundary), 
full regularity for the second derivatives of solutions to the previous systems, 
with Dirichlet boundary conditions (one could also consider slip type boundary 
conditions). The regularity issue for systems like (jl.ip has received substantial 
attention, mostly concerned with an equation in place of a system, and with 
C^°-regularity. In the scalar case, existence and interior integrability of the 
second derivatives are shown in |32| . for any p > 1; in |27j the regularity up 
to the boundary is obtained for any p e (1,2). For systems (solutions are N- 
dimensional vector fields, N > 1), we recall [T] for p e (1,2), [2U] and [33] for 
p > 2, and [22] for any p > 1. These papers deal only with homogeneous systems 
and the techniques, sometimes quite involved, seem not to be directly applicable 
to the non-homogeneous setting. In particular, pQ is the only paper in which 
the L 2 -regularity of second derivatives is considered. However, the results are 
shown only in the interior. Therefore our results seem to be the first regularity 
results, up to the boundary, for the second derivatives of solutions. Another 
main difference with the above papers is that we do not require differentiability 
of S, but merely Lipschitz continuity. For related results and for an extensive 
bibliography we also refer to papers [2J, [H], [T5J, [T7], [IS], [25], [2J], [3D] and 
references therein. 

We have not found papers dealing with the equations arising from the choice 
(|1.4|) for S. This kind of model is used in various branches of mathematical 
physics as, for instance, in non-linear elasticity or in non-linear diffusion. Ac- 
tually, our interest in systems and (|1.3[) arises from our previous studies 
on fluid dynamics problems. Indeed, we recall that a good model for non- 
Newtonian fluids with shear dependent viscosity is the following one 

(1.5) - V- [{[i + \Vu\) p - 2 Vu] +{u- V)u + Vtt= /, V-m = 0, 

which can be obtained from (| 1 . 3[) by adding the contribution of the pressure 
field 7T, the convective term (u- V)it and the divergence free constraint. For this 
system, regularity up to the boundary has been considered in both the cases 
p < 2 and p > 2. The case p = 2 corresponds to the well known Navier- Stokes 
system for Newtonian fluids. For the more general regularity results and a wide 
bibliography on this topic, we refer the reader to [TJ, [PJ for p > 2, and to [SJ for 
p < 2. Despite many contributions to the regularity issue, iy 2,2 -regularity up to 
the boundary for solutions to (ll.5[) is still open, even for the simplified setting of 
"generalized" Stokes system, obtained by dropping the convective term in i|1.5[l . 
We mention the papers [12] and [T3J, which, as far as we know, are the only 
papers where the T^ 2,2 (fi)nC 1,Q (f7)-regularity is obtained, under the additional 
assumption of a small force. The regularity proved below suggests that the main 
obstacle to the iy 2 ' 2 -regularity of solutions of (|1.5I) is actually the presence of 
the pressure term. 

Our interest in fluid-mechanics, and in particular in non-Newtonian fluids, 
leads us to consider the case n — N — 3. However, it is worth noting that our 
results can be immediately extended to dimensions n > 3, and to .ZV-dimensional 
vector fields, N ^ 3. Further, the explicit choices (|1.2p and (| 1 .4[) arc done 
in order to emphasize the core aspects of the results and to avoid additional 
technicalities. Therefore, we do not consider a more general dependence of S on 
Vu or Vu, as for instance S(\7u) = <^(|Vw|) Vu, under suitable assumptions on 
the scalar function tp. For the same reason we avoid the introduction of lower 
order terms. 



2 



In the sequel we cover both the cases p < 2 and p > 2, with, however, some 
differences, and some restrictions on the exponent p , as follows: 

Case p < 2: For p < 2 we consider the "full gradient case" (ll.lj) . In this 
case, all results hold also in the degenerate case /x = 0. For any bounded 
and sufficiently smooth domain f2, we prove W 2,q (£l) regularity, for any q > 2. 
Therefore, we get, as a by product, the Holder continuity, up to the boundary, 
of the gradient of the solution. Results are obtained for p belonging to suitable 
intervals [C, 2), where the constants C are defined precisely. 

Case p > 2: We prove the iy 2,2 -regularity in both cases, and (|1.3D . 

provided that fx > 0. We restrict our proofs to the "cubic domain case" (see the 
next section), where the interesting boundary condition (Dirichlet) is imposed 
on two opposite sides, and periodicity in the other two directions. This choice, 
introduced in reference [5] and used in a series of other papers (see for instance [H 
El El [H] ) > is convenient in order to work with a flat boundary and, at the same 
time, with a bounded domain. The main reason is that, in proving the regularity 
theorem for p > 2 (see Theorem l2.1[) . we apply the difference quotients method: 
we appeal to translations parallel to the flat boundary, and then retrieve the 
normal derivatives from the equations. Then, the simplified framework of a cubic 
domain avoids the need of localization techniques and changes of variables. The 
results can be extended to smooth domains, by following [7], [5], and [§], where 
the extension is done for the more involved system of non-Newtonian fluids (see 
also [2H])- See also the Remark l5Tl 

2 Notation and statement of the main results 

Throughout this paper we denote by fl a bounded three-dimensional domain 
with smooth boundary, which we assume of class C 2 , and we consider the usual 
homogeneous Dirichlet boundary conditions 



Further, we denote by Q the cube Q = (]0, 1[) 3 , and by Y the two opposite 
faces of Q in the a^-direction, i.e. 

r = {x : \x x \ < 1, |ar 2 | < 1, x 3 = 0} U {x : |zi| < 1, \x 2 \ < 1, x 3 = 1 }. 

We impose the Dirichlet boundary conditions on T 



and periodicity, with period equal to 1, in both the x\, x 2 directions. 

By L p (fl) and W m ' p (£l), m nonnegative integer and p G (1, +oo), we denote 
the usual Lebesgue and Sobolev spaces, with the standard norms || • and 
II ' llw m <p(S2)i respectively. We usually denote the above norms by || • || p and 
II ' ||m,p, when the domain is clear. Further, we set || ■ || = || ■ ||2- We denote by 
W 1,p (O) the closure in W 1 *^) of C °°(ft) and by W~ l - p {£1), p' = p/(p-l), 
the strong dual of W ' P (Q) with norm || • ||-i. P '. In notation concerning duality 
pairings, norms and functional spaces, we do not distinguish between scalar and 
vector fields. 



(2.1) 



u\an = 0. 



(2.2) 



u\r = 0, 



We set 



v p {n) = { v G W 



^ p (n) : v ldn = 0}, 
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and 

V P {Q) = { v £ W l p (Q) : u ]r = 0, v is x' - periodic } . 

By Vp(Cl) and Vp(Q) we denote the dual spaces of V p (£l) and V P (Q), respectively. 

We use the summation convention on repeated indexes, except for the index 
s. For any given pair of second order tensors B and C, we write B-C = B^ Cij . 

We denote by the symbols c, c\, c%, etc., positive constants that may depend 
on fi; by capital letters, C, C±, C2, etc., we denote positive constants indepen- 
dent of fi > (eventually, fi bounded from above) . The same symbol c or C 
may denote different constants, even in the same equation. 

We set diU = J^-, <9, 2 u = g ®.g x . ■ Moreover we set (Vu)jj = djUi and 

(T>u)ij — i ((V«)jj + (Vu)ji). We denote by D 2 u the set of all the second 
partial derivatives of u. The symbol D 2 u may denote any second-order partial 
derivative d\ k u except for the derivatives <9f 3 u . Moreover we set 

3 3 

2„, |2 ._ I »2 „. I 2 „„j I r>2„. 2 ._ I c,2 „. I 2 



(2.3) |£ 2 u| 2 := ^l^ttil and | D^u | 2 := ^ |5 



i,j,k=l i,i,k=l 

0,fc)#(3,3) 



We define the tensor 5(A) as 

(2.4) 5(A) = ( M +|A|f- 2 A, 

with [i > fixed constant, p > 1, and A an arbitrary tensor field. It is easily 
seen that 5(A) satisfies the following property: there exists a positive constant 
C\ such that 

(2.5) ^^±B ij B kl >C 1 (v + |A|r 2 |B| 2 , 
for any tensor B. Further 

A - B\ 2 



(2.6) (5(A) - S(B)) • (A - B) > C 2 



(P+IAI + IBI) 2 ^' 
and 

|A-B| 



(2.7) I 5(A) - 5(B) I < C 3 



(M + \A\ + \B\) 



for any pair of tensors A and B, with C2 and C3 positive constants. The proof 
of the above estimates is essentially contained in [21] . We also refer to [16] for 
a detailed proof. 

Our aim is to prove the regularity results up to the boundary given in the 
theorems below. Let us state our main results. We start from the case p > 2. 

Theorem 2.1. Assume that p > 2 and \i > 0. Let f € L 2 (Q), and let u £ 
V^(Q) &e a weafc solution of problem (|1.1[) - (I2.2D or of problem (ll.3[) - (|2.2D . TTien 
it G M /2 ' 2 (Q) . Moreover, there is a constant c such that 

(2.8) ||£> 2 ix||< C ||/||. 
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This theorem will be proved in the next section. 

The other results concern the case p < 2. Note that, in this case, the pa- 
rameter fj, can be equal to zero, thus covering the p-Laplacian systems. Further, 
here we consider a general smooth bounded domain. On the other hand, we 
restrict our considerations to the full gradient case. 

Before stating the regularity theorems for p < 2, let us recall two well known 
inequalities for the Laplace operator. The first, namely 

(2.9) \\D 2 v\\ <C 4 \\Av\\, 

holds for any function v G W 2 - 2 (fl) n W^ 2 {Q.) . Here C 4 = C 4 (fi) . Note that 
if is a convex domain, then C\ = 1. For details we refer to [24] (Chapter 
I, estimate 20). The second kind of estimates which we are going to use for a 
v € w 2 > q {&) n Wo' q (n), q > 2, i 



IS 



(2.10) \\D 2 v\\ q < C 5 \\Av\ 



where the constant C5 depends only on q and f2. It relies on standard estimates 
for solution of the Dirichlct problem for the Poisson equation. Actually, there 
are two constants K\ and K2, independent of q, such that 

(2.11) K iq <C 5 <K 2 q. 

Similarly, one has 

(2-12) || v\\ 2 , q < C||A«||,, 

where the constant C depends on q and tt. For further details we refer to [23] 
and 01]. 

For p < 2 our main results are the following. 

Theorem 2.2. Let be /1 > 0, and 1 < p < 2 such that (2 — p) C4 < 1 , where 
Ci is given by (|2.9p . Let f € Lp+i (fl). Then, the unique weak solution u of 
problem (|1.1[) - (|2.1[) belongs to W 2 ' 2 (il). Moreover, there is a constant C such 
that 



(2-13) IMk 2 <C 11/11 + 11/ 



p+i 



If fi is convex ( or the cubic domain Q ) the result holds for any 1 < p < 2 . 

It is worth noting that in the limit case p = 2, when system (11.11) reduces to 
the Poisson equations, we recover the well known result 

\\U \\ 2 ,2<C 11/11 . 

We set 

(2.14) C 6 =max{C4,C 5 }, 
and 

3g 

(2.15) r(q) = { 3 -(3-g)(2-p) 1 9 < ' 

<7 if g > 3 . 
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Theorem 2.3. Let be fi > 0, q > 2 , and 1 < p < 2 such that (2 —p) C 6 < 1, 
where Cq is given by (|2.14[) . Let / G L r ^ 9 ^(f2) and let u be the unique weak 
solution of problem - (|2.ip . TTien u belongs to W 2,q (Q). Moreover, the 

following estimate holds 

(2-16) \\uh, q <C + 11/11^) . 



Corollary 2.1. Let p, [i and f be as in Theorem \2.3\ . Then, if q > 3, £/ie weak 
solution of problem (|l.ip - (12.1[) belongs to C 1:Q (f2), /or a = 1 — |. 

Note that, in (|2.15[) . r(q) > q for any q < 3. It is worth noting that r(q) tends 
to the same value 3 as q tends to 3 , from below and from above. Furthermore, 
if q = 2 , the estimate (I2.16[) becomes simply (|2.13p . Finally, in estimates (|2.13[) 
and (12.161) . the terms ||/|| and ||/|| q can be replaced by 1. 

Remark 2.1. One could also consider the case where f € L 3 (fl). We omit this 
further case and leave it to the interested reader. Ln this regard we stress that 
our interest mostly concerns the maximal integrability of the second derivatives 
of the solution. 

Remark 2.2. Whenp < 2 we could extend to system (| 1 - 3[) the regularity results 
up to the boundary obtained for system (jl.lj) . by requiring a smallness condition 
on a suitable norm of f . Actually, following arguments already used in \12j and 
for non-Newtonian fluids, the idea is to study the regularity for solutions of 
suitable approximating linear problems and then prove the regularity for solutions 
of the nonlinear problem, by employing the method of successive approximations. 
For brevity, here we avoid this further development. 



3 The iy 2 ' 2 (Q)-regularity: p > 2 and fi > 

In this section we prove Theorem 12.11 Therefore, throughout the section we 
work in the cubic domain Q. Let us introduce the definition of weak solutions 
of both the problems (TLTj) and (TO)) . 

Definition 3.1. Assume that f 6 Vp(Q). We say that u is a weak solution of 
problem (|1.1|) - (|2.2I) . if u E V P (Q) satisfies 

(3.1) / S(Vm) • Vtpdx = I f-ipdx, 
Jq Jq 

for all ip G V P (Q). 

Definition 3.2. Assume that f G Vp'(il). We say that u is a weak solution of 
problem p. 3[) - (12.21) if u G V P (Q) satisfies 

(3.2) / S(Vu)-Vipdx = I f-ipdx, 
Jq Jq 

for all ip G V P (Q). 
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We recall that the existence and uniqueness of a weak solution can be ob- 
tained by appealing to the theory of monotone operators, following J.-L. Lions 

mi- 

In proving Theorem 12.11 we focus on the symmetric gradient case, since the 
full gradient case is, in some respects, easier to handle. Hence we assume that 
S is given by 

S(pu) = {(J.+ \Vu\f~ 2 Vu, 

with fi > and p > 2. 

We follow arguments used in [5J, in the context of non-Newtonian fluids. 
Therefore, we will try to preserve the notations. However in [6] (due to the 
divergence free constraint) the symbol D 2 u has a slightly different meaning from 
that introduced in definition (|2.3[) below, since it also includes the derivatives 
<9f 3 u 3 (see (2.8) in [BJ). 

As in in [BJ , in order to avoid arguments already developed in other papers by 
the authors, we replace the use of difference quotients simply by differentiation. 

It is an easy matter to obtain the following Korn's type inequality, proceed- 
ing, for instance, as in the proof given in |31j . 

Lemma 3.1. There exists a constant C such that 

\\u\\ p + \\Vu\\ p < C\\Vu\\ p , 

for all u £ V P (Q). 

Lemma 3.2. There exists a constant C such that 

||i^u||p<C||V.X>u||p, 

for all u £ V P (Q) . 

This result reproduces Lemma 3.1 in [BJ, adapted to the new definition of 
D 2 u . Note that d s u = on T, s = 1, 2. 

Actually, the above two lemmas hold for each p > 1 . 

Define, for s = 1,2, 

(3.3) J s (u) := / V- [(/i+ \Vu\) p ~ 2 Vu] -d 2 ss udx, 

Jq 1 J 

and 

(3.4) /,(«):= / (M+ \Vu\y- 2 \d s Vu\ 2 dx. 

Jq 

Lemma 3.3. For any smooth function u £ V P (Q) the following inequality holds 
true 

(3.5) J s (u) > d I s (u) , 
with the constant C\ given by (|2.5[) . 
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Proof. Integrating twice by parts in (|3.3p one gets 
J s (u) = f d s \{[i+\Vu\) p - 2 Vu 



■ cLV u dx . 



Note that, due to symmetry, we replace d s Vu by d s Vu. From the above 
expression, one has 

where the derivatives with respect to D^i are evaluated at the point D = T>u. 
Note that here we merely appeal to the chain rule. Then the result follows by 
using estimate (I2.5[) . □ 

Next we prove the following result which, roughly speaking, shows that the 
second tangential derivatives of u are square integrable. 

Lemma 3.4. Assume that f G L 2 (Q) and let u be the solution of problem 
(fT3T) - (j2^|) . Then D 2 u G L 2 {Q) and 

(3.6) P^II<-^rll/H. 

Proof. Multiply both sides of the equations (|l.ip by d 2 s u, s = 1, 2, and integrate 
over Q. By appealing to (|3.3I) and Lemma I3~3l it readily follows that 

/ s (u)<c||/|| \\d 2 s u\\<c\\f\\ ||VM|, 

hence, from Lemma 13. II applied to d s it, 

/.(«)< c || /[| |[0.1>u||. 

Finally, observing that 

^- 2 \\d s Vu\\ 2 <I s (u), 

one gets 

\\d s vu\\< - p _ 2 ||/||. 



Application of Lemma T3. 21 gives the result. □ 

In order to complete the proof of Theorem 12.11 we have to show the integra- 
bility of the remaining second derivatives, namely the normal derivatives <9f 3 u. 
In doing this we follow the argument used in the paper [3]; we express these 
derivatives, pointwisely, in terms of the derivatives of u already estimated, and 
solve the corresponding system in the unknowns 9| 3 Mi, i = 1,2,3. Note that 
the main differences between this situation and that in reference [5], are the 
following: in [3] the L 2 - integrability of <9f 3 u 3 is known thanks to the diver- 
gence free constraint, d 2 3 it 3 = — 9|j ui — d 2 2 1*2- Hence the 3x3 linear system 
considered below is replaced, in [5] , by a 2 x 2 linear system in the unknowns 
d 2 3 Ui, i = 1,2. On the other hand, in reference j3J, the presence of the pressure 
prevents the full IV /2 ' 2 -regularity. 

For the missing derivatives we prove the following lemma. 
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Lemma 3.5. Let u be the solution of problem (|1.3p — (|2.2[) . Then the vector field 
d 33 u satisfies the pointwise estimate 

(3.7) | 8 2 33 u |<c(i p _ a |/|+ | Dlu |) , a.e. in Q . 



Proof. Straightforward calculations show that 

d s [(p+ \Vu\) p - 2 Vu] = (p+ \Vu\f- 2 d s Vu 
l "" S ' +(p-2)0 + \Vu\) p - 3 \Vu\- 1 (Vu- d s Vu)Vu. 

For convenience, we set Vj k = (Vu)j k and B := (fx + \Du\). By using (|3.8D . 

the j th equation for any j = 1, 2, 3, takes the following form 

(3.9) 

BP- 2 (d 2 kUj + 8 2 k u k ) + - 2)BP- 3 \Vu\- 1 V lm V jk (d 2 kmUl + d 2 kl u m ) = -2f r 



Let us write the previous three equations as a system in the unknowns d 3 
For j — 1 , 2 we have 



33 u j 



(3.10) 5 P ~ 2 4s «i + 20 - 2) 5 p ~ 3 \V u]- 1 V j3 ]T V l3 8 2 33 ui = Fj —2 f 3 , 

i=i 

where 

2 3 



F,:= -BP- 2 J2d 2 kkUj -BP- 2 J2d 2 k u k 



k=l k=l 

(3.H) 3 

-2(p-2)BP- 3 \Vu\- 1 d 2 kmUl V jk V lm . 

l,m,k = l 
(m,fc)^(3,3) 

For j = 3 we have 

3 

(3.12) 2£P- 2 9 3 2 3 Uj + 20 - 2) BP- 3 \Vu\~ 1 V j3 ]T P /3 d 2 33 u t = F j - 2 

i=i 

where, for j = 3, 

F 3 -= -BP- 2 j^d 2 kk u 3 -BP- 2 j^d 2 k u k 

k=l k=l 

(3.13) 3 
-2(p-2)BP- 3 \Vu\- 1 £ dl m u x V ih V lm . 



l,m,k — l 



The equations (|3.10|) . for j = 1,2, together with the equation (13. 13)) for j = 3 
can be treated as a 3 x 3 linear system in the unknowns d 33 uj, j — 1,2,3. 
Multiply all three equations by B 2 ~ p . We denote the elements of the matrix 
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A = A(x) associated with this system as a,ji, where j, 1 — 1,2, 3. Then, we can 
write the system in a compact form as 

(3.14) a 3 i <9f 3 ui = Gj , 

where the elements of the matrix of the system are given by 
a jt :=5 fl + 2(p-2) (B \Vu\y 1 V j3 V i3 , 

for j = 1,2, by 

a jt := 2Sji + 2(p - 2) (B \Vu\y 1 V ]3 V l3 , 

and for j ' = 3, and 

(3.15) Gj := B 2 ~ p (Fj - 2 fj ) . 

Note that ciji — aij] moreover, if £ denotes any vector field then 



= | £ | 2 + £f + 2 (P - 2) (B \Vu\ y 1 [Vu ■ £ 



Hence, the matrix A — (<Xji) is also definite positive, a.e. in x 6 Q, and the 
previous identity shows that 

«,/ i, o > lei 2 - 

By setting £ = 9| 3 u, we have obtained 

(3.16) \di 3 u\ 2 < \G\\di 3 u\, a.e. in Q , 

where, obviously, by G we mean the vector (G\,G2,G 3 ). Noting that, from 
(|3~T5)1 . (j3~TT|) and (IXT3)) . there holds 



(3.17) \ Gj l< Tp _ 2 \ fj \ + c \Dl u l a .e. in 



2 

from this estimate and (|3.16[) we get (|3 . 5|) . □ 
Finally, by combining (|3.6p and Q3.5p we readily obtain 

P 2 "ll<^- 2 II /II 
which is just (|2.8[) . The proof of Theorem 12. II is accomplished. 

4 A regularity result for an approximating sys- 
tem: p < 2 . 

In the sequel we introduce an auxiliary positive parameter r/ and study the 
regularity for solutions of the following approximating problem 



(4.1) 



-r/Ai)-V'5(Vu) = /, inn, 
v = 0, on dfl , 
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with S defined by (|2.4|) . r\ > 0, \i > and p £ (1, 2). The solutions satisfy the 
estimate (14.151) below, with the constant C independent of r\. This allows us to 
show that, as r\ —> 0, v v tends, in a suitable sense, to the solution v of problem 
(|4.1j) with rj = 0. A similar situation occurs, with respect to fj,, as /i — >• 0. 

We explicitly note that we introduce the above model just to approximate 
our solution by smooth functions. 

Let us introduce the definition of weak solution of both the problems (|4.1[) 
and (fl~il-(l2~TD. 

Definition 4.1. Assume that f £ V^'(f2). We say £/ia£ tj is a weak solution o/ 
problem (|4.ip i/u € a7id satisfies 

(4.2) 77/ Vu-V<^ote + / 5(V t>) • dx = / f-(pdx, 

Jn Jn Jn 

for all ip £ V2 (fi) - 

Definition 4.2. Assume that f £ K,'(f2). We say that u is a weak solution of 
problem p.ip - (|2.ip . if u £ V p (fl) satisfies 

(4.3) / S(Vu) • V</>cte = / f-ipdx, 
Jn Jn 

for all ip £ V p (il). 

As recalled in the previous section, the existence and uniqueness of a weak 
solution is known from the theory of monotone operators. 

We start by proving the W 2,2 -regularity result stated in Proposition 14.11 
below. In (|4.4j) . the dependence of the constant c on Qo, 77 and /j is omitted 
since the aim of the proposition is just to ensure that second derivatives are well 
defined a.e. in Q. Following the notations introduced in section 2, by capital 
letters, C, C\, C%, etc., we denote positive constants independent of /j and of 77 
also. 

Proposition 4.1. Let p £ (1,2), / £ L 2 (fl), and v be a weak solution of 
problem (|4.1I) . Then v £ W,' (CI) and, for any fixed open set Hq CC fl, there 
exists a constant c such that 

(4-4) || D 2 v\\ LHQo) <c 11/11 . 

Proof. As in the previous section, we formally use derivatives instead of differ- 
ence quotients, to make the computation simpler. Fix an open set Oq CC O. 
Let £ be a Cq (O)-function, such that < C( x ) < 1 in 51, and C( x ) = 1 in fio- 
Multiplying the first three equations in (|4.1j) by — V • (£ 2 V v) and integrating 
over n we get 
(4.5) 

77 I d?jVid h (( 2 d h Vi) dx+ I dj [{n+ \\7v\) p - 2 (Vv) l3 ] d h (( 2 d h Vi) dx 
Jn Jn 

= - fi d h (C 2 d h Vi) dx . 
Jn 
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By integration by parts, with respect to Xj and Xh, on the left-hand side one 

has 

(4.6) 

V f (d? h Vi) 2 { 2 dx+ f d h \Wv\f- 2 (Vv)^] d h (Vv) l3 ( 2 dx 



n 



= -r) / {d 2 jh Vi) Ri jh (x)dx- / d h [(p+ \Vv\)p~ 2 (Vu)ij] R ijh (x)dx 
Jn Jn 

- I U{d 2 h v l )C 2 dx-2 f f i (d h v i )C(d H Qdx=Y,I i , 
Jn Jn i=1 

where, with obvious notation, Rijh are lower order terms satisfying estimates 
(4.7) \Riih{x)\ < c|C||VC||V»|. 

As in the proof of Lemma T3. 3 1 it is easy to verify, by appealing to (|2.5[) , that 
(4.8) 

[ d h \Vv\)P- 2 (Vv) l3 ] d h (Vv) l3 C 2 dx> c [ (p+ \Vv\) p - 2 \D 2 v\ 2 C 2 dx. 

Jn Jn 

On the other hand, by Holder's and Cauchy- Schwartz inequalities, 

(4-9) |/i|< e |!I^ICII 2 + C ( £ )||vci! 2 00 liv«|| 2 , 

(4-10) |/ a |< e\\\D 2 v\t\\ 2 +c(e)\\f\\ 2 , 

and 

(4.H) \h\< c||VC||oo||/||||Vt;||. 

Further, by using the estimate 



< c(^+ \A\y- 



dA kl 
we have 

(4.12) |/ 2 |<c f(fx+ \\7v\) p - 2 \D 2 v\ |C| |V C||Vw| dx , 

Jn 

and, by the Cauchy-Schwartz inequality, 
(4-13) \h\< 6|||^|CH 2 + C (e)||VC||Ll|V«|ir 

From (|4.6p together with ||Vw|| < c||/||, it follows that 
(4.14) II \D 2 v\ C||<c H/ll . 

Hence fl~4l holds. □ 

Our next step is to get a global estimate for the L 2 -norm of the second 
derivatives, uniform in rj. This is the aim of the following proposition. 
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Proposition 4.2. Let be (2 - p)C 4 < 1, with C4 given by (|2.9I) . Let / £ 

L5+r(fi), and /et v be a weak solution of problem (|4.ip . TVien w belongs to 
W 2 ' 2 (fl). Moreover, there exists a constant C such that 



(4-15) IMk2<C II/H + 11/11 



p+i 



Proof. In order to avoid a useless dependence on /z, we assume, without loss of 
generality, /1 £ (0, 1]. At first note that, by replacing ip by v in (14. 2[) it is easy 
to get the following estimate for ||Vu|| p , uniformly in 77, 



Since, by Proposition 14.11 v E ' (f2) 5 the « th equation (|4.1I) can be written 
almost everywhere in f2 as 

r?A^+ ( M + |Vw|) p - 2 A^ 

+(p - 2)(/i + |Vv|) p - 3 |V»| _1 V« • (dj V«) 9j m = -fi- 

By multiplying both sides by Auj and summing over i = 1, 2, 3, we have 

n\Av\ 2 + |Vv|) p " 2 |Av| 2 

= (2 - + |Vv|) p-3 |Vu| -1 Vu • (dj Vv) dj Vi Av, - f t Av % , a.e. in Q . 

Next, we drop the term ?7|Au| 2 , and bound the left-hand side from below by 
(fjt+ \Vv\ ) p ~ 2 |Au| 2 . Multiplying the estimate thus obtained by (fi+ |Vv|) 2_p 
and then integrating over Q we get 

f I Av \ 2 dx < (2-p) / I L> 2 v I I Aw I da; + / (p + |V«| ) 2 ~ p | / 1 1 Av \ dx , 
Jn Jn Jn 

where we have used the estimate (for details see the Appendix) 

I Vv ■ {dj Vv) (dj v t )A Vl \< \Vv\ 2 \D 2 v\\Av\. 

Observing that (jjt + | Vv|) 2 ~ p < fi 2 p + |Vw| 2 p , using Holder's inequality, and 
dividing both sides by || Av ||, we get 

(4.17) ||A«|| <(2-p)\\D 2 v\\ + \\\S7v\ 2 -rf\\+ || / || . 

Let us estimate the first two terms on the right-hand side. For the first term 
we employ estimate (|2.9[) . As far as the second term in (|4.17[) is concerned, by 
applying Holder's inequality with exponents 3/(2— p) and 3/(p+l), the Sobolev 
embedding of W 2 ' 2 (fl) in W 1,6 (£l), and by appealing to the estimate (|2.12[) with 
q = 2, we get 

II |V V | 2 - p /|| < || S7V || 2 - p || /|| a < C II A V || 2 ~ p || /|| a . 

P+1 P+1 

By using the above estimates in (14.170 , we get 

II Av || < (2 - p) C4 Av\\+C || Av || 2 - p || /|| a + || / || . 

P + 1 
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Recalling that (2 — p)C± < 1, and applying the Young's inequality 
(4.18) a 2 - p b<ea+ c(e)b^ , 

it is easy to recognize that the estimate 



(4-19) II A« || <C (II/H + ll/H V 



p+i 



holds. By using once again (|2.12p we prove (14.151) . □ 

5 The iy 2 ' 2 -regularity result: p < 2 . 



Proof of Theorem \2. 6 A We deal separately with the case fi > and the degen- 
erate case (i = 0. 

The case ji > - Consider the "sequence" (v^) consisting of the solutions 
to problem (|4.1|> . for r\ > 0. By the above proposition the sequence (v v ) is 
uniformly bounded in W 2 ' 2 (fl). Therefore, by Rellich's theorem, there exists a 
field u € W 2 ' 2 (Q) and a subsequence, which we continue to denote by (v v ), such 
that Vjj — ^ u weakly in W 2,2 (£l), and strongly in W 1,q (fl) for any q < 6. Let us 
prove that 



(5.1) 



/ S( Vu) • V<p dx = lim < / S(Vv„) ■ Vipdx + r\ \ Vv v -Vipdx> , 
Jn v^o+ [J Q Jn ) 

for any ip G C^°(f2). By applying (12.71) and then Holder's inequality, we get 

/ S( Vu) • Vip dx - I S(Vv v ) - Vipdx 
Jn Jn 

<c (n+ \Vu\ + \Vv v \ ) p ~ 2 I Vu - Vv v | \V<p\ dx 
Jn 

<c [ | Vu - Vv v | p_1 \V<p\ dx <c\\ Vv v - Vu W*- 1 || Vtp \\ p . 
Jn 

The right-hand side of the last inequality tends to zero, as r\ goes to zero, thanks 
to the strong convergence of v v to u in W 1,p (il). Further 



/ Vv n ■ V<p dx 
Jn 



<77||Vi;„|| ||V^||, 



where the right-hand side tends to zero as 7/ goes to zero. Finally, observing 

that for any 77 > and any <p e Cq°(D,) the right-hand side of (|5.1I) is equal to 

J n f-ipdx, we show that u satisfies the integral identity (|4.3p for any ip 6 Cq°(CI). 

By a standard argument we show that u satisfies the integral equation (j4.3p , for 

any ip 6 V P (Q). Hence u is a weak solution of (jl.ip . and belongs to W 2,2 {£l). 

Moreover, (12 . 13[) follows from the relation || u |[a 2 < liminf ||u^ || 2,2 , together 

77— >o+ 

with (I4.15p . From the uniqueness of weak solutions we obtain the desired result. 

TTie case fx = - Let us denote by u M the sequence of solutions of (II. ip 
for the different values of /1 > 0. We have shown that the sequence (u M ) is 
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uniformly bounded in W 2 ' 2 (£l). Therefore, exactly as above, we can prove 
the weak convergence of a suitable subsequence in W 2 ' 2 (fl), and the strong 
convergence in W l q (£l) for any q < 6, to the solution u £ W 2,2 (£l) of the 
problem with /i = 0. In this regard note that estimate (|2.7|) also holds 

with [i = 0. 

Finally we prove the last assertion in Theorem 12.21 For a smooth convex 
domain ft estimate (|2.9[) holds with C4 = 1. Hence the assumption on p is 
merely p > 1 . □ 

Remark 5.1. We could adapt the above arguments to the case p > 2. Via 
a result similar to Proposition \4-l\ one shows that the solution v of the ap- 
proximated system \4--l\ belongs to Wj ' c (f2). Then reasoning as in the proof of 
Proposition \4-.2\ one obtains a global estimate for v in W 2 ' 2 (fl), uniformly in r\, 
with a restriction on the range of p, p £ (2,2 + ^-), C4 as in (|2.9p . Hence, as 
Theorem \2.2\ above, one proves that the solution of with (i > 0, belongs to 

W 2 - 2 (fl). This result has the advantage to be directly proved in a general smooth 
domain, without need of localization techniques. However, it requires limitations 
on the range of p and, moreover, it cannot directly cover the case fj, — 0, since 
the W 2 ' 2 (fl)- estimates that one obtain are not uniform in /1. 



6 The VF 2 ' 9 -regularity result: q > 2 and p < 2 . 

Proof of Theorem \2.3l From Theorem 12.21 we already know that the solution u 
of problem (jl.ip belongs to W 2 ' 2 (il), since (2 — p)C^ < 1 . Therefore, we can 
write equation (|4.16p with u in place of v, and rj = 0. By multiplying this 
equation by (/i + |Vw|) 2 p , we can write, a.e. in ft, 

i « . , Vu ■ VVu • Vw „, ._ ,.0-d 

^ - Au -^ Vr^y^ = /( ^ +|Vu|) ' 

where we have used the notation Vu ■ VVw • Vu to denote the vector whose i th 
component is Vu ■ (dj Vu) dj Ui = (di it/.) (d^ Uk) (dj Ui) . 

We start by proving an a priori L 9 -estimate for the second derivatives of u by 
assuming, for the moment, that u £ W 2,q (Vl). We follow an argument similar 
to that used for proving the V^ 2,2 -estimates of u. We multiply both sides of 
equation (|6.ip by — Au |Au| 9 ~ 2 , and integrate in Q. We get (for details see the 
Appendix) 

f \Au\ q dx< (2-p) [ \D 2 u\\&u\ q ~ 1 dx+ f {ii+\Vu\) 2 -v\f\\&u\ q ~ 1 dx. 
Jn Jn Jq 

By appealing to Holder's inequality and to the inequality (fi + |Vu|) 2 ~ p < 
1 + |Vu| 2 ~p, we show that 

||A«||« < (2-p)||C 2 u|| 9 ||Au||^ 1 

+ ll/ll 9 l|Au||«- 1 + |||Vu| 2 ^/|| 9 ||Au||«- 1 . 

Further, by dividing both sides by || Au one gets 

(6.3) || Au \\ q < (2 - p) || D 2 u ||, + || / 1|, + || |Vu| 2 ~ p / ||,. 
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We estimate the first term on the right-hand side of (I6.3[) via inequality (|2.10p . 

Concerning the last term on the right-hand side, we start by assuming that 
q G (2, 3). As usual we denote by q* — 3p/(3 — p) the Sobolev embedding 
exponent of q. By applying Holder's inequality, with exponents s = q*/(2 —p)q 
and s' — r(q)/q, we get 



(6-4) |||V U | 2 ^/|| 9 <||V W ||^||/| 



r{q) 



From (|6.3p . by appealing to (12.121) , (|6.4|) and to Young's inequality one easily 
gets 

[[Au||, < (2 -p) C 5 \\Au\\ q + ||/|| 9 + e ||Au||, + c(e)||/||^ . 
Recalling the assumption on p, a further application of estimate (|2.12[) gives 



(6.5) \\u\\ 2 , q <c[\\f\\ q + ll/n;- 

Next we assume that q > 3. We will use arguments similar to the previous 
ones. Actually, by appealing to the Sobolev embedding W 1,q (£l) L°°(tl), to 
the estimate (|2.12p , and to Young's inequality, we estimate the last term on the 
right-hand side of (|6.3p as follows: 
(6.6) 

II \Vu\ 2 - p f\\ q < || fj> II /||, < C || An || 2 ^ || /||, < e \\ Au\\ q + c(e) \\ f\\f" . 

Then, by repeating verbatim the arguments used above, one shows that u is 
bounded in W 2 ' q (fl), uniformly with respect to /i, and that the estimate (I6.5P 
holds. Finally, the argument used in the proof of the Theorem 12.21 in order to 
extend the results to the degenerate case \i = apply here as well. 

The previous arguments are formal, since we have assumed that solutions 
belongs to W 2 ' q {VL). However the following argument applies. Let us consider 
the problem 
(6.7) 

w £ = , on <9fi , 

where w £ is the unknown and J £ denotes the Friedrichs mollifier. The coefficients 
of this modified system belong to C°°(lR n ). We can also write this system in 
divergence form, as follows: 

(6.8) -d h [m ijhk (x)d k w e ]+( P -2)d h [c e ijhk (x)] 3 k w £ = f ( M + |V U |) 2 ~ P , 
where 

m ijhk (x) = SijShk + {p - 2) c\ jhk {x) 

and 

c ijhk( x ) = dhJe{u % ) d k J e { Uj ) 

J (A* + Js(\ VW|)) Je(|Vu|) 

Further, let 

c ijhk (x) = (d h u i )(d k u j ) {fM+ ^ ul) |Vtt| , 
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From the well known estimate 



(6.9) |VJ e (u)| = |J e (Vu)| < J e (|Vw|) 
we get 

(6.10) \ c ljhk( x )\ — 1) uniformly in x , £ , and \i . 

This shows that the system (|6.8[) (hence, the system (|6.7I) ) is a linear elliptic 
system with regular coefhcients. For such a system it is well known that if a force 
term F belongs to L q (il), q > 2, then the solution belongs to W 2,q {Vi) (see, 
for instance, [H]). By following the previous arguments with u replaced by w e , 
and by using (|6.9I) and (|6.10p , it is straightforward to obtain the estimate (|6.5[) 
for w e . Note that such estimates are uniform with respect to \x and e. Hence, 
there exists a subsequence, still denoted by w £ , and an element w £ W 2 ' q (tt) 
such that, as e goes to zero, w e converges to w, weakly in W 2,q . Convergence 
is also strong in W 1,r (0,): for any r if q > 3, and for any r S ^1, if q < 3. 
Let us show that w is a solution of the system 

. , . . , Vw • VVw • Vw . , ._ ,,0-d 

^ - Aw -^- 2 ) (, + |V,|) |Vu| =/ ^ + ' Vu ') • 

To this purpose, we write equations (|6.7|) and (16. lip in the weak form, and take 
their difference, side by side. This leads to the expression 

(d h wf - d h w t ) d h (fi dx + (2 - p) I (c e ijhk - c ijhk ) d\ k w] tp l dx 
^ 3.12) ^ Jn 

+ (2 - p) / Cijhk (dhk w j ~ d lk w j) ¥>* dx > 
Jn 

for any cp S C^°(ft). The first integral goes to zero as e goes to zero, thanks to 
the strong convergence of w £ to w in W 1,2 (il). Concerning the second integral, 
we recall that mollificrs converge in L p to the mollified function, as e goes 
to zero, and that LP convergence implies almost everywhere convergence of a 
subsequence. Therefore, c\^ hk converges to Cijhk ; a - e - m ^- From (I6.10p . by 
recalling that fi is bounded and by using the dominated convergence theorem, 
it follows that 

(6.13) lim / \cf jhk - c ijhk \ 2 dx = . 

e ^°Jn 

Hence the second integral in (I6.12j) goes to zero. The last integral in (|6.12[) 
tends to zero, thanks to the weak convergence of w e to w in W 2 ' q (fl), since the 
coefficients c^hk (x) are bounded. 

Finally, it is easy to verify that w = u . Indeed, by taking the difference of 
(|6.ip and (|6.1ip . side by side, and by setting V = u — w, we get 

= 0, in Q , 




Finally, multiply the above equation by AV and integrate in Q. By appealing 
to arguments already used, one readily recognizes that, under our assumptions 
on p, the vector V satisfies || AV || = 0. Hence V = 0, by uniqueness. 

□ 
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The Corollary 12. II is an immediate consequence of Theorem 12.31 Details are 
left to the reader. 



7 Appendix 

Our aim is to show the estimate 

1 1 := | V« • (dj Vv) {dj Vi) An, | < | Vv| 2 \D 2 v\ \Av\ . 

In the sequel, for convenience, we sometimes avoid the summation convention, 
by explicitly writing the sums, even if repeated indexes appear. 
We recall that 



3 3 

and \D\\ 2 :=<T(D 2 Vk) 2 '= E \ 9 >* 

k,j,h=l 



j,h=l fe=l 

We introduce the vectors b and w, whose components are defined as follows 

3 

bj := (dj v) ■ Av , w 2 k := E ( ( d h «fc) h 3? ■ 
The modulus of vector b satisfies the following estimate: 

3 3 3 3 

h = E b > < = \ Av \ 2 ED^) J = i a -i 2 i^i 2 . 

j = l j=l j = l i=l 



Hence 
(7.1) 

Moreover 



3 3 

'jt = J2(d h v k ) 2 J2b^\Vv k \ 2 \Av\ 2 \Vt 

3=1 



h=l 



3,h,k=l 



3 

^E 

fc=i 



E (dhVk)bj 
3,h=l 



3 

^E 



E 

=i \ j,fc=i 



* E (( d hVk)bjf, 



where, in the last step, we have used that, for any pair of tensors A and B, 
there holds \A ■ B\ < \A\ \B\. Hence, by the above notations and estimate (|7.ip . 
we get 



i / 1 < E i^-i ^ \ Av \ \ Vv \ E ' Vwfc i 



fc=i 



< \Av\ IVwI 



fe=i 



\ Ei^i 
\ fe=i 



\ El V ^l 2= |At>| |V«| 2 |D a v| 
\ fe=i 
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which is our thesis. 
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